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We demonstrate entanglement génération between mode pairs of a quantum field in a single, rigid 
cavity that moves nonuniformly in Minkowski space-time. The effect is sensitive to the initial state, 
the choice of the mode pair and bosonic versus fermionic statistics, and it can be stronger by orders 
of magnitude than the entanglement dégradation between an inertial cavity and a nonuniformly 
moving cavity. Detailed results are given for massless scalar and spinor fields in (1 + 1) dimensions. 
By the équivalence principle, the results model entanglement génération by gravitational effects. 
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In the emcrging ficld of relativistic quantum informa- 
tion, the past décade has revealed a number of novel kine- 
matic effects on entanglement between observers in rela- 
tive inertial motion [ï| and between one inertial and one 
uniformly linearly accelerating observer (see Refs. [IHsl 
for a sélection and Réf. Q for a récent review). The 
phenomenon that has received most attention is the en- 
tanglement dégradation caused by relative accélération, 
assuming that the issues of switching the accélération on 
and off arc negligiblc and ficld modes of plane wave type, 
sprcad over large régions of the space-time, are consid- 
ered. A rare case of entanglement génération under thèse 
assumptions is found in Réf. [ÏÔj . 

A framework that removes the idéalisation of spatially 
unlocalised field modes is to confine the field to a cav- 
ity [lH, and the idéalisation of everlasting accélération 
can be removed by taking the cavity's worldtube to bc 
inertial outsidc a finite time interval [Ï2j |. This créâtes 
a situation that should be more directly cxpcrimentally 
accessible. When the cavity's worldtube consists of iner- 
tial and uniformly linearly accelerated segments, it was 
shown in Refs. 12, [ï]| that the évolution of the quan- 
tum field is amenable to a systematic perturbative anal- 
ysis when the accélérations are small compared with the 
cavity mode frequencies, both for bosonic and fermionic 
fields; further, the perturbative analysis remains valid for 
arbitrary durations of the individual segments, and the 
distances travellcd may hence be arbitrarily large. The 
nonuniformly accelerated motion was found to dégrade 
mode entanglement between the moving cavity and an 
inertial référence cavity. The physical mechanism bchind 
the dégradation is that when observations of the field are 
restricted to a small number of field modes, for exam- 
ple by appealing to the frequency sensitivity profile of a 
detector by which the field might be observed, particle 
création within the moving cavity causes information to 
be lost into the infinitely many field modes that are not 
bcing observed. 

In this article wc show, within the above perturbative 
treatment of a quantum ficld that is confincd to a single 



cavity, that nonuniform motion of a rigid cavity créâtes 
entanglement between any two given field modes within 
the cavity, even when the initial state is separable. We 
develop a gênerai quantitative analysis for both a scalar 
field and a fermion field, giving detailed results for a sam- 
ple travel scénario and demonstrating that the particle 
statistics has a significant effect on the entanglement. 
As a highlight, wc show that the entanglement généra- 
tion can appear in the first order in the small accéléra- 
tion expansion, while the entanglement dé gra dation be- 



tween cavities is only a second-order effect [12J, ll3j . This 
suggcsts that entanglement génération in a cavity may 
be more readily observable than the entanglement dégra- 
dation between cavities, allowing more effective tests of 
phenomena linked to the dynamical Casimir effect 



and possibly also improving the technological prospects 
of building quantum gâtes based on accélération effects. 
For instance, highly entangled two-mode squeezed states, 
produced by a known gâte in continuous variable Systems, 
can be generated by periodically rcpeating segments of 
uniformly accelerated and inertial motion [15|. 

We work in (1 + l)-dimcnsional Minkowski space with 
metric signature ( — h): additional transverse dimensions 
can be included via their contribution to the effective 
field mass. The length of the cavity in its instantaneous 
rest frame is ô > 0. The cavity is assumed to be inertial 
outsidc a finite time interval, but the initial and final 
vclocities need not coincide. We use units in which h = 
c = 1. Complcx conjugation is denoted by an asterisk 
and Hcrmitian conjugation by a dagger. 0(x) dénotes a 
quantity for which 0{x)/x is bounded as x — > 0. 

Bosons: Wc consider a real scalar field (fi satisfying the 
Klein-Gordon équation (— □ + m 2 )ifi = 0, where m > 
is the mass and □ is the scalar D'Alambcrtian. For def- 
initeness, we adopt at the cavity boundaries the Dirich- 
let boundary condition as in [12j, although much of the 
analysis holds for any boundary condition that ensurcs 
unitarity of the time évolution. 

Let {<fi n | n = 1, 2, . . .} bc a complète orthonormal set 
of mode solutions that are of positive frequency with rc- 



2 



spect to the cavity's proper time at early times (the in- 
rcgion), and let {</>„ | n = 1,2,...} be a similar set at late 
times (the out- région) . Each set has an associated set of 
laddcr opcrators, with the respective commutation rela- 
tions [a^a^] = 5 nm and [à„,à]„] = 5 nm , and a vacuum 
state, denoted respectively by | ) and | } . The two sets 
of modes are related by the Bogoliubov transformation 



and the ladder opcrators are related by 
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where the notation is as in Réf. [ldl. The vacua are re- 



lated by (17 1 
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where W := \ J2 pq Vpc 
normalisation constant 

Wc prépare the System in the in-rcgion in a statc with- 
out mode entanglement. We ask: does the cavity's mo- 
tion gcnerate mode entanglement in the out-region, where 
the particle content of the state has changed? 

Our mcthodology is as follows. We first specify the 
in-region state and express it in the out-region basis us- 
ing ([3]) and the adjoint of (J2j) . We then trace over ail 
out-region modes except those labelled by two distinct 
quantum numbers k and k'. We quantify the entangle- 
ment in the resulting reduced density matrix by the neg- 
ativity [Ï8l42Ô| , defined as minus the sum of the négative 
eigenvalues of the partial transpose. The advantages of 
negativity are that it is easy to compute and it interpo- 
lâtes betwcen entanglement monotones that have a more 
direct operational interprétation 21 1 . 

To idcntify a parameter régime that is treatable ana- 
lytically, we assume that the Bogoliubov coefficients have 
Maclaurin expansions in a small dimcnsionless parame- 
ter h, such that 



a = a<® +a« + a< 2) + 0(h 3 ) 
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where the superscripts indicate the power of h, a 1 -- ^ = 
diag(G*i, G2, • • • ) and each Gj has unit magnitude. For 
cavity worldtubes grafted from inertial and uniformly ac- 
cclcratcd segments, this situation arises when the accél- 
ération at the centre of the cavity is proportional to h/ S 
by a numerical coefficient that may differ from segment 
to segment: the parameter h is in this case the product 
of the cavity's width S and the accélération at the centre 
of the cavity [I2I ], We then work perturbatively in h. It 



follows that to order h 2 we have N = 1 — j J2 P q | 2 
and 

|o> = (i-ŒKT)|ô) + Œ^X|ô) 

pq pq 
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As a first example, we take the in-region state to be the 
in- vacuum |0). To order h 2 , the partially transposed, 
reduced density matrix vanishes outside a 6 x 6 block. 
Among the six eigenvalues, the only possibly négative 
ones are 
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where /£_,„ := \ Y. q ^ n \Pqra 
der h 2 . A4 arises from cohérence between 



and Vfcfc' is kept to or- 
Ô) and 

| lfc ) | ïfe'), while A6 arises from cohérence between | Ô) 
and \ %)\2 k ,). 

Spccialising to a cavity worldtube that is grafted from 
inertial and uniformly-acccleratcd segments, we find that 
a qualitative différence émerges depending on the relative 
parity of k and kl . If k and k' have opposite parity, the 
expansions given in Réf. [Ï2I ] and their massive generali- 
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in the negativity is then linear in h and given by \i~> kk i\- 
If, by contrast, k and k 1 have the same parity, we have 

$i = and Vk^ = v^) + °( h3 )- The leadin g term 

in the negativity cornes then from A6 and is of order h 2 . 
Sample negativity plots for both cases are shown in Fig. 
[T] for a massless field when the cavity undergoes a single 
segment of uniform accélération. 

As a second example, we take the in-rcgion state to 
be |lfe), containing cxactly one in-particle. Using 
and the adjoint of ((2]), we find 
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To order h 2 , the partially transposed, reduced density 
matrix now vanishes outside an 8 x 8 block. Among the 
first five eigenvalues, the only possibly négative one is 



j kk ' 
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which arises from cohérence between | ïfc ) and 
| 3fe ) | 2^/ ). The last three eigenvalues are the roots of 
a cubic polynomial, analytically cumbersomc for generic 
values of the parametcrs but readily amenable to numer- 
ical work. 

Specialising to a cavity worldtube that is grafted from 
inertial and uniformly-accelerated segments, we again 
find a qualitative différence depending on the relative par- 
ity of k and kl . In particular, if k and k' have opposite 
parity, the lcading contribution to negativity cornes from 
the eigenvalue 



*-kk' 
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and is linear in h. The negativity is in this case higher 
than the corresponding negativity for the in-region state 
| 0). Sample negativity plots are shown in Fig. [T] for a 
massless held when the cavity undergoes a single segment 
of uniform accélération. 

Fermions: We consider in the cavity a Dirac field ip 
with mass m > 0, with boundary conditions that ensure 
unitarity of the time évolution as in Réf. [l3j . 

Let {ip n | n e Z} be a complète orthonormal set of 
mode solutions, such that the modes n > are of posi- 
tive frequency and the modes with n < are of négative 
frequency with respect to the cavity's proper time in the 
in-region. We write the in-region mode expansion of the 
field as tp = J2n>o ^rJon + J2n<a ^ci, so that the non- 
vanishing anticommutators of the fermionic ladder opera- 
tors for positive frequency modes and négative frequency 
modes are {& m ,6t} = 5 mn and {c m ,c) n ) = <5, 
tivcly and the in-vacuum || )) satisfies b n || 
c« I )) — 0. The corresponding modes, operators and 
states in the out-region arc denoted by a tilde. 

We write the Bogoliubov transformation between the 
two sets of modes as ip m — J2 m A m n VVu where the Bo- 
goliubov coefficients A mn form a unitary matrix. The 
relations between the two sets of ladder operators can be 
written in terms of the Bogoliubov coefficients by tak- 
ing appropriate inner products. The vacua are related 
by || » = Me w ||Ô)), where W = £ P >o V M SJ4« M is 
a normalisation constant, and V pq can be expressed in 
terms of A mn [Ï3| . 

We assume again that the Bogoliubov coefficients have 
a Maclaurin expansion, now of the form 



mn respec 
and 
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where h is the small dimcnsionlcss parameter, A 
G m ô mn (no sum) and each Gj has unit magnitude. We 
then work perturbatively in h fl3j ]. 

We first take the in-region state to be the in-vacuum 
I )) and trace over ail out-region modes except a mode 
k > 0, of positive charge, and a mode k' < 0, of néga- 
tive charge. As pointed out in Réf. [22}, the tracing in a 




(b) 

FIG. 1. The leading order contribution to the negativity is 
shown for a massless scalar field and a massless Dirac field. 
The travel scénario has a single accelerated segment, of accél- 
ération h/ S as measured at the cavity's centre and of duration 
t in the cavity's proper time, where the corresponding Bo- 
gogliubov coefficients can be found in Refs. P3, • The neg- 
ativity is periodic in u := hr/[A8 atanh(/i/2)] with period 1. 
Fig. Q^a) shows M /h, in dashed (blue) for a scalar field with 
in-region vacuum and (k,k') — (1,4), in dotted (blue) for 
a scalar field with in-region state | lfe ) and (k, k') = (1,4), 
in solid (red) for a Dirac field with in-region vacuum and 
(k, k!) = (2, —1) with s = 0, and in dotted-dashed (red) for a 
Dirac field with in-region state || 1 K )) and (k, k) = (1,4) with 
s — in the notation of [TJ- Fig. QJb) shows the correspond- 
ing curves for N/h 2 with the scalar field modes (k, k') = (1, 3) 
and the fermionic modes (k, k!) = (1, —1) and (k, k) — (1, 3). 



fermionic Fock space has an ambiguity, but the construc- 
tion of V guarantees that this ambiguity does not affect 
our measures of entanglement, given that thèse measures 
are unambiguous for the in-region state || )) . To order h 2 , 
we find that the only potentially négative eigenvalue of 
the partially transposcd, reduced density matrix is 
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f A , + f A , 
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Note the similarity 



where 
with (|6b|) . 

Specialising to the massless fermion and a cavity 
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worldtube that is grafted from inertial and uniformly- 
accelerated segments [Ï3| , we find that the eigenvaluc (fTTT) 
is linear in h whcn k and k' have opposite parity. The 
leading order correction to the negativity is then equal to 

\A%\. 

Consider then in-region states with particles. When 
the in-region state contains a single particle in mode 
k > 0, the out-region reduced density matrix for two 
modes of opposite charge turns out to have vanishing 
negativity to order h 2 . However, when considering two 
modes of the same charge as the excitation, e.g., k, k > 0, 
cntanglcmcnt is gencrated between thèse modes. Mirror- 
ing thc vacuum case an analogous expression to Eq. (jlip . 
with f£-, K i and f^ K replaced by and f£, is ob- 

tained for the possibly négative eigenvalue of the partially 
transposed, reduced density matrix. Similar as before the 
leading order correction to the negativity is | A^l | if k and 
k have opposite parity. 

This structure is simply a conséquence of the fermionic 
algebra and charge conservation: entanglcment is gen- 
erated by the Bogoliubov coefficients Connecting two 
modes, either by création of a pair of oppositely charged 
particles in thèse modes, or by shifting a préexistent exci- 
tation to another mode of the same charge. If the mode in 
question is populated, the Pauli exclusion principle pro- 
hibits any further excitation of this mode and thc création 
of particle-antiparticle pairs cannot entanglc it with any 
mode of opposite charge. 

Whcn the in-region state contains a pair of oppositely 
charged particles, the situation changes again: the par- 
tial transpose of the out-region reduced density matrix of 
the corresponding modes k > and k' < then has one 
potentially négative eigenvalue of the form of Eq. (fTTj) . 
where f^ K > an d f^-, K a re replaced by f£ and f^,. Spc- 
cialising to the massless fermion and a cavity worldtube 
that is grafted from inertial and uniformly-acccleratcd 
segments [Ï3|, we find that when k and k' have oppo- 
site parity, the leading order correction to the negativity 
is the same as when the in-region state is the in-region 
vacuum. 

Samplc plots are shown in Fig. [TJ 

Conclusions: We have demonstrated that nonuni- 
form motion of a cavity générâtes entanglement between 
modes of a quantum field confined to the cavity, both 
bosonic and fermionic. Working to quadratic order in 
the cavity's accélération, and quantifying the entangle- 
ment by the negativity, we found that the entanglement 
génération dépends on the initial state of the field, on the 
relative parity of the mode pair that is observed at late 
times, and on the bosonic versus fermionic statistics. For 
both bosons and fermions, we found situations where the 
entanglement génération can be enhanced by placing par- 
ticles in the initial state. For fermions, however, charge 
conservation and the Pauli exclusion principle require the 



choice of the considered out-region modes to be consistent 
with thc initial state to generate entanglcment, while thc 
bosonic statistics allow the modes to be freely populated 
without hindering cntanglcmcnt génération. 

Comparcd with the motion-induced entanglement 
dégradation between a static cavity and a moving cav- 
ity 12,01, wc found that the entanglement génération 



can occur already in linear order in thc cavity's accélér- 
ation, while the entanglement dégradation is a second- 
order effect. The prospects of expérimental vérification 
[ÏH could hence be significantly better for phenomena 
signalling entanglement génération than entanglement 
dégradation. 

Thc motion-induced entanglement effects that we have 
analysed have technical similarities with the création 
of squeezed states in resonators with oscillating walls, 
known as the dynamical Casimir effect [H, 23|. In this 
context, wc emphasise that our only approximation was 
to work in the small accélération régime, meaning that 
the product of thc cavity's length and accélération is 
small compared with the speed of light squared [TU, [Hl ■ 
Our analysis hence covers as a spécial case cavities that 
oscillate rapidly with a small amplitude: such cavities are 
often introduced in theoretical analyses of the dynamical 
Casimir effect but are experimentally problematic (Ï3 |. 

Our analysis however covers also cavities that accel- 
erate in a given direction for finite but arbitrarily long 
times, with travel distances that may be arbitrarily large. 
Further, as the équivalence principle implies that gravi- 
tational accélération can be locally modcllcd by accéléra- 
tion in Minkowski space-time, our results suggest that 
a gravitational field can producc entanglcment. Ex- 
periments for cntanglcmcnt génération could hence be 
sought in setups that span macroscopic distances, includ- 
ing quantum communication through near-Earth satellite 
orbits. 
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